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O 1 I. INTRODUCTION 

<N 

Oh' The multiscattering problem is hard to solve in a fully quantum-mechanical context. Consider, for instance, inclusive 
pion-nucleus scattering, which will be our model system. Around the A(3, 3) resonance several channels are open, 
besides the elastic one: absorption of the pion, inelastic, single charge exchange and double charge exchange. Some 
reactions can take place several times for the pion inside the nucleus. All the reactions interfere with each other 
and their typical reaction probabilities are strongly dependent on the region of the nucleus. Cascade methods 0] 
which reduce the complicated looking output of the reaction to simple steps, seem appropriate to deal with such a 
problem. However, cascade methods often involve drastic semi-classical simplifications not completely under control. 
\ Typically the reaction probabilities are taken from the free cross section or at most Pauli blocking is included (and 
f^S ■ as a consequence are one-body mechanisms) and a classical propagation is used in between two collisions. Even so, 
\f~} ' the potential ability to embody known physics of the problem and the versatility of the cascade approach makes it 
, interesting to study how it could be improved to systematically include quantum corrections. 

Cascade method have been widely used in the context of heavy ion collisions, where one wants to describe the 
dynamical evolution of a large number of particles. In the crudest approach only two-body collisions of free particles 
are considered, without any mean field effect ||. Complementary to this is the hydrodynamic approach or the 
more sophisticated time-dependent Hartree Fock method, where the mean field is well described but collisions are 
neglected g[. Both models are merged into microscopic kinetic or transport models. Particularly successful has 
been the Boltzmann-Uehling-Uhlenbeck (BUU) equation which embodies mean-field effects, two-body collisions and 
Pauli blocking [1 ulpj. Transport mode ls g o back to Boltzmann and are widely used in physics of plasmas of all 
kinds ||. (See [L0[ pd|^[T^, [l3| [l6| [17], [l^] for reviews on the transport theory approach from different points 

• rH , of view.) In heavy ion collisions, microscopic transport models are applied not only to nucleons but also to mesons 
(in particular pions) and resonances [fl9|| . Transport models have also been used directly for studying pion-nucleus 
H \ reactions |2l], |2^, |24|, |2^, |2(| or other meson-nucleus reactions |27j . 

In classical kinetic theory the basic quantity is the distribution function in phase-space, f(x,p,t), which satisfies a 
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transport equation of the type gain- loss as in Eq. (1.3) below. It has long been recognized that the Wigner transform 
pM p9| is the natural way to derive a quantum transport equation || [n], ||(], EH]]. The subject has been brought 
to a high degree of sophistication, as needed for instance in the description of quark-gluon plasmas, where it has 
to include thermal effects, be consistent with relativistic invariance as well as covariance under non-Abelian gauge 
transformations [l4[ ^3[ [m| |35| . The field of quantum transport theory is currently quite active and with 
many open issues, both of conceptual and of practical interest 

The mark of cascade models is the classical propagation of the particles between two collisions. This, which at first 
sight may be considered as a drawback, could also be one of their main virtues, for this quasi-classical propagation 
does indeed exist, furthermore it is dominant between to well separated successive collisions. Consider, for instance, 
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the following form of the uncertainty principle 

|At| ^ \F F ~ Pn^T' (Lla) 

|Ax| Pa , — T . (1.1b) 

' ' IP-Pon-ahell^)! 

These relations show that only on-shell states can propagate a long time or distance, and that these long lasting states 
can be taken as classical in the sense that only they survive as fi if formally taken to zero. Furthermore, from Eqs. 



(1.1) one finds for large At, Ax 

A-r rJW, 

(1.2) 
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on-shcll 



which is the classical Hamilton's equation. This crude argument can be made more precise once a definite prescription 
is chosen to simultaneously use conjugate variabl es s uch as (t, E) or [x,p) in the quantum mechanical context, for 
instance, Wigner's prescription. Of course, Eqs. ( |i~l| ) only hold as long as interaction does not take place, but they 
show the existence of two scales in the multiscattering problem: one due to the mean free path of classical states 
and the other due to the uncertainty principle, i.e., propagation of virtual states. Multiscale problems usually make 
trouble to approaches which fail to explicitly include such a feature. It seems more promising to divide the problem 
into two parts: first, the quasi-local virtual states are integrated out in such a way that the elementary vertices of 
the microscopic theory are substituted by classical effective ./V-body quasielastic plus absorption probabilities, and 
second, these probabilities are then used in a cascade approach, where only quasi-classical states show up explicitly. 

In order to clarify the meaning of the rather abstract program above, we should place it in the appropriate context. 
One of our motivations has been to understand and justify the success of the approach in |24|, [25]. There all the 
pion- nucleus inclusive reactions are computed along the following lines: the pion self-energy is computed in nuclear 
matter including all the Feynman graphs considered relevant for energies around the resonance. The imaginary part 
of this self-energy is then considered as a pion "width" against decay of the elastic channel into reaction states. 
With the help of a local density prescription, the total reaction cross section is then calculated. Furthermore, use is 
made of Cutkosky rules ^| in order to separate the reaction width into the several reaction channels, namely, 
absorption and quasielastic (with or without charge exchange). This information is then used in a Monte Carlo 
simulation of the path of the pion inside the nucleus: the pion is treated classically in between collisions but the 
reaction probabilities are computed microscopically from Feynman graphs. Several questions arise from this rather 
intuitive approach, such as how does it follow from a purely quantum-mechanical calculation?, how can the classical 
and local density approximations be systematically improved? how to avoid double counting between, say, a genuine 
three-body absorption mechanism included in the microscopic calculation and a quasielastic followed by a two-body 
absorption coming from the Monte Carlo simulation, having both the same final state? The Monte Carlo simulation 
for the pion can be cast in the form of a transport equation, 

d t f(x,p,t) = J d 3 x' d 3 p> Q(x,p;x' ',p')f(x' ,p' \t) - R(x,p)f(x,p,t) , (1.3) 

where Q(x,p; x',p') represents the unit time probability for the transition (x',p') — > (x,p) in phase space and R(x,p) 
is the probability of leaving the state (x,p). Q contains both elastic and quasielastic processes while R is the total 
reaction rate. The absorption rate is thus 

A{x, p) = R{x, p)- d 3 x' d 3 p' Q{x', p'; x, p) . (1.4) 

Then, which are exactly the kernels R, Q, if any, that will produce the same cross sections as the Schrodinger 
equation does? Note that we shall actually deal with just one particle cascading through the nucleus so f(x,p,t) 
is the probability density in phase-space rather than the real density as in standard kinetic theory of plasmas. The 
equation is however formally identical to a transport equation. 

In this work we address those issues. Our main concern has been to write exact quantum-mechanical equations 
in such a way that the connection with cascade methods were immediate. To this end, a density matrix formalism 
is used plus Wigner's prescription, in order to achieve a well-defined classical limit. The Wigner transformation has 
been applied before to study scattering, most notably by Remler in a scries of papers |58| |j| flfj, |l| . Our emphasis 
is, however, different since we are interested in isolating the ./V-body quasielastic and absorption rates seen by a 
quasi-classical particle so that the correct quantum results are recovered in a cascade model. 
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We shall not be concerned here with relativistic invariance (although relativistic kinematics is allowed and actually 
used in applications involving pions at intermediate energy), thermal effects, or gauge invariance, however, we shall 
find that some of the findings in these more sophisticated fields are also of interest here. We shall find that a full 
Wigncr transform, in space and also in time, is needed in order to include inelastic channels, a point not usually 
realized in non-relativistic applications of the Wigner transform to collision theory. 

As we have said practical cascade models are often very simplified from the many-body point of view. Typically, one 
deals with classical particles moving in a mean field potential, and classical collisions using in-vacuum cross sections 
or decaying width in-vacuum lifetimes, restricted by Pauli blocking in the final states. This procedure mimics at a 
classical level the evolution described by Feynman graphs. The problem with a direct diagrammatic approach is, of 
course, that a realistic description would require to carry out the computation to graphs of arbitrarily high orders. 
The cascade method aims at an efficient procedure to carry out a resummation of those graphs, at the price of a 
classical description. To improve on this approach it would be interesting to make a formulation using exact relations 
between resummed sets of graphs, much in the line of the kinetic Kadanoff-Baym or the Schwinger-Dyson equations, 
but with an explicit h dependence and within a space-time framework, and then find a prescription to integrate out 
the virtual, non-classical, intermediate states. The improvement of this approach resides in the fact that, once precise 
definitions for the effective reaction probabilities are given, one is no longer constrained to use in-vacuum estimations 
for them and, in principle, many-body effects can consistently and systematically be included. This work should be 
regarded as an attempt in this direction. Within a particular simplified model we try to implement the previous 
program. In doing so, we introduce ideas, some of them hopefully new, that presumably will be present in future, 
more systematic, developments. 

The paper is organized as follows: In Section || we show that indeed knowledge of the evolution of the density 
matrix in the Wigner's f orm , taken as a density of classical projectiles, gives rise to the correct quantum-mechanical 



cross section. In Section III the one particle system is studied with emphasis on its Wigner's form and the classical 
limit. In Section |y| the same analysis is carried out for a particle in an optical potential, i.e., for the elastic channel. 
Section [v] is devoted to the general many-body evolution equation exemplified with a simple pion-nucleus model, with 
nly pions and particle-hole (ph) excitations as physical degrees of freedom. Section |vi| sho ws how the ph degrees of 



o 



freedom can be removed in order to obtain a purely pionic evolution equation. In Section VII the virtual pionic de grees 



of freedom are identified and integrated out in order to obtain effective TV-body quasielastic rates. In Section VIII 



the actual outcome of our scheme is illustrated in simple cases. Finally in Section IX we summarize our conclusions 



II. CROSS SECTION IN THE SIMULATION APPROACH 

In this Section we shall assume that the proper simulation procedure has already been carried out (how to do that 
will be the subject of subsequent Sections) and our present purpose is to show that in this case the correct fully 
quantum mechanical cross section is obtained by the usual method. We start by relating the S'-matrix in Wigner's 
representation to the cross section, and later contact will be taken with the time evolution operator which is closer 
to an actual simulation procedure. 

Let us assume that the Monte Carlo simulation gives us directly the pi n , p ou t relationship 

Pout = S 5' t , (2.1) 
where p is the density matrix and S is the scattering matrix. Then we shall have 



Pout(u) =/ d b vS{u,v)p in (v). (2.2) 

Here u, v are points in the phase space (x,p), d 6 v = d 3 x v d 3 p v , p(u) is the density matrix in Wigner's form |^8| , |29| , 
00 

p(x,p) = J d 3 y e-*y-P' h (x + \y \p\ x - \y) (2.3) 

and S(u,v) is a real function related to the Wigner's form of S to be interpreted below as the "probability" density 
of going from v to u in phase space. 

The cross section from some initial state \i) to some final state |/) is (see for instance p2l) 

<7(i-/)= / d 2 b {f\Se- iP - b l h \i)\ (2.4) 
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b being the impact parameter vector and P the momentum operator. In addition \i), /} represent normalized initial 
and final states [}42[ . Using the property of the Wigner's representation 



ti(XY) 



(ft 



■X(u)Y(u) 



the cross section can be rewritten as 



tr(»->/) = J d 2 bd 3 xd 3 p d ^ P 3 p f (x',p') S(x',p';x,p)pi(x - b,p) 



(2nh) 

Now for I/} we take a plane wave with momentum pj 

pf(x,p) = (2irhf 5{p - p f) . 
This, of course, spoils the dimensional counting, but that will be fixed later: 



a(i — > pj) — J d 2 bd 3 x' d 3 x d 3 p S(x' ,p f ;x + b,p)pi(x,p) 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



The next step is to use that the fine details of t he p rojectile wave function are not relevant, i.e. only its momentum 
distribution matters. This allows us to rewrite (2.8) in the form 



a(i—>p f ) = J d 2 bd 3 x d 3 pS(x' ,p f ;b,p) J d 3 xp l (x,p) . 



Equation (2.9) follows from Eq. (2^8) provided 

V x J d 2 bd 3 x S(x',p f ;x + b,p) = 0. 



(2.9) 



(2.10) 



Physically this is clearly true: the integration on b projects out the x± dependence. Furthermore, varying scii amounts 
to a change in the initial position of the projectile along the same incoming trajectory, but that only shifts x' which 
is integrated out. A proof of (2.10) is given in Appendix |A]. 

Then for \i) normalized but with a narrow momentum distribution we can take 



and so 



J d 3 xpi(x,p) = pi(p) w (2irh) 3 8{p - p^ 

a( Pl -> p f ) = (2nh) 3 J d 2 bd 3 x' S(x',p f ; b, Pi ) . 
The correct dimensions are recovered by integrating p * within a solid angle dtt / 

! >o(Pi^p f ) = J dpfp 2 d 2 bd 3 x / S ( x ^ pf . bp j 



dClf 



(2.11) 



(2.12) 



(2.13) 



This is the desired relationship between cross section and Monte Carlo output: to obtain da/dQ. several projectiles with 
random impact parameter, represented by the integration over 6, should be thrown against the target, the simulation 
procedure (contained in the function S(x' ,Pf]b, Pi )) will put them in (x',Pf) after the interaction, but only the 
scattering angle is relevant (represented by the integration over x' and pf). S(x' ,Pf,b, Pi ) can then be identified with 
the density probability of going from (b,Pj) to (x',Pf) due to the interaction. It can be noted, however, that this 
"probability density" is not necessarily positive. Positivity is only required for the cross section, i.e., after integration. 
Likewise in general, conservation of energy p J = \ Pi \ is only achieved after integration over b and x' . 

Another remark is in order regarding ( 2.13 ). Unitarity of S (or also Eq. ( |2.2[ )) implies 



d 6 u S(u, v) 



(2.14) 
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then a direct use of ( pTl3| ) yields 



r 



dfi 



f 



d b — oo . 



(2.15) 



As usual, this means that the non-interacting part of S must be removed from the cross section: 

S(u,v)=Z(u)8(u-v)+T(u,v), (2.16) 
where T is less singular than S(u — v) and represents the scattering probability while Z(u) is the non-interaction 



probability. T instead of S should be used in ( 2.13 ) (whether S is unitary or not). 

The S matrix is a convenient theoretical tool but actually the simulation procedure is more directly related to the 
evolution operator rather than the S matrix: a projectile is sent and the outgoing distribution is observed much later. 
This is described by 



Pl2 



e -i(t 2 -H)H/h ^ i(t 3 -ti)H/h 



-h,t 2 ^+oo. (2.17) 
In Wigner's form we then have 

p(u, ta) = J d 6 v U(u, v; t 2 ~ ti)p(y, h) . (2.18) 
It has to be shown that for large t = t-x — ti, U(u, v; t) does the same job as S(u, v) did in ( 2.1 3| ) . For large enough t 

e —i(t 2 -ti)H/h ^ e ~it 2 Ho/h g e itiH /h _ e -itH /h g ^ 

where Ho, H represent the free and full Hamiltonians respectively, and we have used the well known property 
[Hq, S] — [^2|, This allows to relate the two functions S(u, v) and U(u, v; t) as in Eq. (A. 2), namely, 

U(x' )P ';x,p;t)n I f^e^'-y^-'iMP'HQhMP'-H^/ns^p'-^p). (2.20) 
J (2nhy 



Now, if U is used in (2.13) instead of S, the x' integration gives q = and the result is the same, that is, 



f 



J dp fP 2 f d 2 bd 3 x'U{x',p f ;b,p i ;t). 



(2.21) 



Once again, the non -inte raction probability is to be removed from U. Also note t hat th e result does not depend on t 
since q vanishes in ( 2.20 ) after x' integration. This requires t to be large so that ( 2.19 ) holds. 

Let us summarize the outcome of this Section. A simple-minded method to obtain the cross section would be: 
1) to compute the kernel of the evolution equation, U(u,v;t); 2) to use it in a simulation procedure as a transition 
probability density from the point v in phase space at time t\ to the point u at £2; 3) to count the outgoing particles 
to extract the cross section as in a real experiment. What has been shown here is that the naive expec tation is indeed 
correct. Note that we use "simulation procedure" to mean a method solving an equation like ( 2.18 ). The fact that 
U(u, v; t) will not be positive definite in general can be a technical problem but the equation itself is well-defined. 



III. THE SINGLE-PARTICLE EVOLUTION EQUATION 

In this Section we revi ew t he simplest case of a one-particle system and its connection with a classical description. 
Let us construct Eq. ( |l.3D for a single-particle system, with Hamiltonian H. We start with the density matrix 
evolution equation 



iti 



dp 
'~dt 



H,p 



(3.1) 



and rewrite it using the Wigner transformation, as defined in (2.3), to eventually consider its classical limit. The 
product of two operators can be dealt with by means of the identity 



(AB) (x,p) 



d 3 yd 3 q d 3 zd 3 k 
(nh) 3 {nh) 3 



A(y, q) B(z, fc) e «(9- fc W« e »(M>)-W/ft e «(p-«l)-*/s . 



(3.2) 
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where the left-hand side stands for the Wigner's form of AB at (x,p). This can be written more compactly using the 
notation (simplectic scalar product) 



u A v 



-v A u 



xq-yp, u=(x,p), v = (y,q) 



as 



(AB)(u) 



(ttTi) 3 (ttH) 3 



(3.3) 



(3.4) 



In order to study the classical limit, we should transform this expression into one with better properties as fi goes to 
zero. To this end we use the identity (in one dimension) 



e rx P /h = 2irhe lRa * a "6{x)5( P ) 



(3.5) 



where the exponential in the right-hand side is to be expanded as a series of powers of h. This identity can be 
established by considering both sides as distributions on e lax as test function, for arbitrary a. It can immediately be 
extended to any number of dimensions, and also 



iuAv/h 



6 ihd^Kdv 



S(u)S{v) . 



(3.6) 



Then (3T) can be cast in the form 
(AB)(u) = 

or simply p3| 

Expanding in powers of fi yields 



A(u + v)B(u + w)(TTh) 6 ei lhd ^ d ™5(v)8(w) 



(irh) 3 (nh) 3 
ei iha » Aa <"A(u + v)B(u + w) 



v—w—0 



(AB)(u) 



ih 



(AB) («) - A(u)B(u) + - {A(u), B(u)} P + ■ ■ 



(3.7) 



(3.8) 



(3.9) 



The zeroth order shows that operators commute in the classical limit and the first correction introduces the usual 
Poisson bracket, {A, B} p — dA A dB. 



The evolution equation for p(u,t), Eq. (3.1), takes the form: 



ih 



9p(u,t) 

at 



e 2 



H(u)p(u,t) 



At lowest order in Ti we find 



2ism(\hd {H ^ H(u)p(u,t) . 

dp 



Of 



{H,p} p + 0{h z ), 



(3.10a) 
(3.10b) 

(3.11) 



which is the classical equation of evolution in phase space. This equation is also referred to as the Liouville equation in 
mechanics, the Vlasov equation in plasma physics (particularly when the particles move coupled to an el ectrom agnetic 
field), or (collisionless) transport equation in the context of kinetic theory. Correspondingly, Eq. ( 3.10b ) is the 
quantum transport equation. The approach of using the Wigner transformation to derive a quantum transport 
equation has become standard [ pT| and has been extended in various ways, including finite temperature [ fli] , |32| | 
relativistic treatments fl6), Abelian and non-Abelian gauge covariant definitions of the Wigner function as well 
as second quantization definitions of a Wigner operator []10j| . 

For subsequent developments, it will be convenient to consider also the more general case of a non-Hermitian 
Hamiltonian H — Hr + iHi. The new evolution equation is 



ihf t =Hp-pm. 



(3.12) 
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and the corresponding equation in Wigner's form would be 

d tP (u, t) = | cos (\fid [H) A d (p) )H I {u)p{u, t) + H sin (±?id (H) A d ip) )H R (u)p{u, t) 

= ^Hi{u)p(u,t) + {H R {u),p(u,t)} P + 0(n). (3.13) 

Throughout this Section we have used the classical limit in a rather formal manner. The problem of fixing the h 
dependence in an expression is not a trivial or even well-defined one. In practice, our point of view has been to adopt 
the Wigner's form of an operator as that with a smooth classical limit or plainly as that without any Ti dependence 
at all. However, it is clear that this cannot be true for all operators |62|. After all h is not small, it is rather unity 
in natural units. The classical limit should be understood as a physical limit of small fluctuations. Depending of 
the system this may correspond to large times or distances, large momenta or energies, weak coupling (or sometimes 
strong coupling), low densities or large number of degrees of freedom, among others. 



IV. EVOLUTION EQUATION AND OPTICAL POTENTIAL 



In this Section, we extend the study of the evolution equation for the density matrix and its classical limit by 
including some many-body effects by means of an optical potential. Let us consider the following energy-dependent 
Hamiltonian: 

H(E) = H + V op t(E), (4.1a) 

J 7T tj — bj' + IT] 

Hq is energy independent and Hcrmitian, and contains the free and Hartree-Fock pieces of the particle self-energy. 
On the other hand, V op t(E) is energy dependent and hence non Hcrmitian and non instantaneous, and contains the 
intermediate states not in the elastic channel. iVi denotes the anti-Hermitian part of V op t- This is the absorptive part 
of the optical potential and is non positive definite. 
The Schrodinger equation in this case takes the form 

E\E) = H(E)\E) , (4.2) 

or in time representation 

ihd t \if>,t) = J dTh(T)\if),t-T), (4.3) 

where 

NM> = / ^e-' B,/ V(B)l-B) , (4.4) 

- S(t)H q + ^-9{t) I dE e- iBr / n Vi(E) ■ (4.5) 

7Tft ./ 



The non- locality in time in (O) comes from non-elastic intermediate states represented by the optical potential. As 
in the previous Section, we would like to construct an evolution equation for p(t) = \ij},t)(ij},t\, however, it is easy to 
see that p(t) does not satisfy an autonomous equation [|63| . The problem is that p(t) does not actually contain the 
same information as t) does; it loses track of the phases. This would be irrelevant if \tp, t) where the wave function 
of the whole system, but it is not: it only describes the elastic part. Instead, we have to consider the more general 
set of operators 

p(t,t') = \ip,t)(ip,t'\ (4.6) 
as well as their Wigner's transformed (in time-energy) 

p(t, E) = J dr e iET ' h p{t +\r,t- ±t) (4.7) 
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which are Hermitian. Note that the equal-time density matrix is recovered through energy integration: 

f dE , 



(4.8) 



Similarly we define the Wigner transform of the Hamiltonian H (with H(t,t') := h(t — t')) 

H(t,E)=H(E). (4.9) 
In order to write an evolution equation for p(t,E), it is more convenient to start with the energy representation: 



E\E) = H(E)\E) 
p(E u E 2 ) = |£i)<£ 2 | 



(4.10a) 
(4.10b) 



so that 

EkftEuEa) =H(E x )p(R L ,E i ), (4.11a) 
E 2P {E 1 ,E 2 ) = p(E u E 2 )HHE 2 ). (4.11b) 

Subtracting both equations and using the notation E — \{E\ + E 2 ), co = E\ — E 2 and p(u, E) — p(E\, E 2 ), we obtain 



up(uj, E)=H{E+ ±co)p(u, E) - p(oj, E)H^ (E - |w) , 
which can also be written in terms of p(t, E) by Fourier transforming u>: 

ihd t p(t,E) = e? lha t P>d z H) H{E)p{t,E) -e-i^ 9 ^ p(t,E)H Jf {E). 



(4.12) 



(4.13) 



This is the evolution or transport equation in presence of an optical potential. It is an autonomous equation for p(t, E) 
as a function of t because, as a consequence of energy conservation, no energy derivatives of p appear in it. Note that 
if we ta ke th e "classical" limit in the right-hand side by taking the explicit % — > 0, Eq. (3.12) is recovered. 

Eq. (4.13) has been written using a partial Wigner form, namely, in time-energy. It can be further expanded by 
using the full space-time Wigner's transformation (in time-energy and po sition- momentum) . This yields the following 
more symmetrical form of the transport equation which generalizes Eq. (3.10a) 



ihd t p(u) = e-^ H)Ad ^ P) H(u)p(u) - e^ H)Ad ^ p(u)H^u) 



(4.14) 



where now u denotes (t, x; E,p), and u A v = tE' — Et' — x • p' + p ■ x'. This formula holds actually for a general 
non conservative and non instantaneous Hamiltonian H(t,E). In what follows we will not, in general, expand the 
formulas by expliciting the position- momentum part of the Wigner transform. 

As it stands, ( 4.13 ) is of little usefulness since it is non-local in time and so all the time derivatives of p will 
contribute to the first o ne. Fo rtunately, the expansion in powers of h can be used to bypass this undesirable feature. 
First let us expand Eq. ( 4.13 ) using the separation H(E) — H + Vr(E) + iVi{E), 



ihdtp = 



Ho + V R ,p 



9eV r , y<9 t/ S 



iVi, p 



+ 



id E Vi, 



ih\ 2 



(4.15) 



where [,]_,[, ] + stand for commutator and anticommutator respectively. Th en, le t us demand that there exist a 
classical limit at all for dtp(x,p; t, E). The first term in the right-hand side of ( 4.15 ) poses no problem because it is 
of order Ti due to the commutator (cf. ( |3.9[ )). Analogously the third and higher order terms have explicit Ti in them. 
So we will require 



Vi = o(h) 



(4.16) 



for dtp to exist in the limit U — * 0. Note that this implies Vr = 0{K) too, due to ( 4.1b ). Under this assumption, it 
follows that the first and the second terms in Eq. (4.15) are 0(h). The third term is 0(fi ). The fourth term is 0(h 3 ), 
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and the others are 0(h 4 ) and 0(h 3 ). In general, higher order time derivatives of p appear only at higher order in h. 
This allows to express d t p in terms of p (t, E) only (without time derivatives) at any given order in Ti: d t p in the third 
and fourth ter ms ar e eliminated using (4.15) recursively, dfp in the fifth and sixth terms are eliminated by applying 
{ifidt) once to (4.15), etc. Explicitly, through second order, we obtain 



ifidtp = H ,p + i Vi,p 



+ 





1 












V R ,p 


- + 2 


d E V R , 


H ,p 


+ i 


Vi,p 


+_ 



(4.17) 



This expansion is expected to be only asymptotic. In what follows we will not, in general, explicitly expand the 
equations to put them in a manifestly instantaneous form, as done here, but this procedure can be carried out if 
needed. 

It is worth noticing that the real part of the optical potential Vr, only appears at higher order than the free 
Hami ltonian Ho or the absorptive part of Vj. In many cascade calculations only Vj is used, (i.e., the cross section). 
(4.17) indicates that this is a kind of classical approximation. 

As we have seen, the relation (4.16) is needed in order for Eq. ( [4. 17] ) to make sense. Its origin is clearer if V opt is 
considered as a self-energy in many-body language [j47| . From Eq. ( |4.l| ) the Hartree pieces are included in Hq and 
then V op t contains only self-energy pieces with loops, the only diagr ams w ith imaginary part. As shown in |4q| , each 
loop gives a further power in H to a diagram, consistently with Eq. ( 4.16 ). 

Let us consider the propagation of a particle in infinite nuclear matter. Due to translational invariance, all relevant 
operators ar e func tions of the momentum and commute among them (internal degrees of freedom are neglected here) . 
In this case ( [4.17 ) becomes 



d t p(t,E) = 



T(E) 



p(t,E) 



(4.18) 



where T(E ) — 0(K) is the width and is given as a power series in h. A closed form is more easily obtained directly 
from (|T|) 



-iT(E) = H(E - ±iT) -H\E+ \iT) . 
From here, expanding in n, the well-known quasiparticle result Wq 



T(E) a -2Z(E)Vi{E), 
Z(E) = (I - d E V B {E)y l 



(4.19) 

(4.20a) 
(4.20b) 



appears in order fi . 

It can be noted that from the two equ ation s (4.11) we have extracted only one equation, namely, (4.12). Taking 
the semi-sum instead of the difference in ( 4.11 ) one obtains 



Ep{t,E) 



_ e ^hd t d E H(E)p(t,E) + -e"^ im 9 



- p{t,E)W{E). 



(4.21) 



which is a kind of energy-shell constraint equation |], [w). At lowest order in fi, in the full Wigner form, it just says 
that E = H(u ), i.e., the density matrix is on-shell. By construction this equation is consistent with the transport 
equation ( |4.13 ). In fact, using again an asymptotic expansion in fi, we can put this equation in instantaneous form. 
The lowest orders are 



Ep 



Ho,p 





i 














- + 4 


d E Vi, 


H ,P 


+ i 
















+_ 



+ o(h 3 



(4.22) 



If the transport and const raint equations are written as dtp = L t p and Ep — LeP, with the linear actions L t and 
Le defined by ( 4.17 ) and ( 4. 22] ), respectively, the compatibility amounts to the statement [L t , Le] = 0, which can be 
verified also by explicit calculation to the order shown. 
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As often emphasized, the two equations, transport and constraint, are needed for a proper description of the 
evolution || |l(| |l4|, |3^, |34|. In next sections we will also make use of two equa tions, however, our approach will 
involve a differ ent set of kinematic equations, ( [5.13 ) and ( 5.15 ), or ( |6.11 ) and ( 6.1 3| ) ) . This is further discussed at the 
end of Section |Vl| 

An equal-time density matrix is often used in non-relativistic transport equations |38 . In the relativistic case the 
equal-time formulation can also be used 
transform As discussed in 1 33 



but Lorentz invariance is only manifest by using the space-time Wigner 
both formulations are equivalent without introducing further approximations, 
the space-time formulation being nevertheless richer since it contains the energy distribution |33|, |3~4| . 

The equivalence of the two approaches (equal-time vs. time-energy Wigner distribution functions) holds whenever 
the Hamiltonian is instantaneous, that is, whenever H (u) is independent of E. This follows from integrating over the 
energy in (4.14) (a procedure named energy averaging in p3L p4j), recalling the relation (4.8) between the two density 



matrices, and noting that 8e acts only on p(u) (it is not necessary to assume that the Hamiltonian is conservative, 
H(u) may depend on i). The situation for a non instantaneous dynamics is different and it requires the use of the 
space-time form, as done here. This is not related to relativity but rather to the fact that the particle can leave the 
clastic channel and spend some time in other inelastic intermediate states before returning to the elastic Hilbert space. 
In this case to derive an autonomous evolution equatio n for the equal-time density matrix is no longer straightforward. 



We will deal with a closely related problem in Section VII 



V. NON-ELASTIC STATES 



In the previous Section we studied the time evolution of the density matrix describing the elastic part of the wave 
function. Here we would like to study how to describe the non-elastic part. To fix ideas consider the scattering of 
pions by nuclei in the A-isobar region |5(| . Typical processes are those depicted in Fig. [l]. Fig. [j]a shows a self-energy 
graph contributing to the pion optical potential: the incoming pion tt collides with a nucleon of the nucleus which is 
excited to a A-isobar and leaves a hole h in the Fermi sea of nucleons, the A couples further with a Nir' state. If 
these intermediate states are only virtual, i.e., for a short time of order Ti, the pion n" emerges with the same energy 
as 7r: it is an elastic scattering and it only contributes to the real part of the optical potential. On the other hand, if 
the intermediate particles ir'Nh are near their mass-shell, the \tt' Nh) state can live a long time and we have instead 
the process in Fig. [l]fr: it is a real decay of a pionic mode into pion-particle-hole and contributes to the imaginary part 
of the optical, i.e., gives a width to the "elastic" incoming pion tt. 

We wish to describe such quasielastic reactions (Fig. [!]&) by means of a density matrix formalism appropriate to 
connect with cascade calculations. In order not to unnecessarily obscure the discussion we shall use a simplified model 
with two kinds of particles: "pions" and "ph" (particle-hole), both bosons, without explicit A-isobar or isospin degrees 
of freedom. In addition, the equations here will not include pion absorption. The absorption mechanism is developed 
in Appendix [S] and added later. The model is given by the following Hamiltonian 

H = + tf (ph) + Hj := H f + Hi (5.1) 
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a) 



b) 



FIG. 2: The same processes as in Fig. pi, for the model in Eq. (5.1) 



where H$\ H^ h) are one-body operators for pions and ph, and Hi is the interaction vertex 7T7rph: 

Hi = J d 3 xd 3 yd 3 zF{x;y,z)ij) ph (x)4>l(y)^(z) + h.c. 
:= F + F\ 



(5.2) 



and as usual 



V) 0ph 



= o. 



(5.3) 



The interaction vertex is assumed to be elementary (i.e., instantaneous and without energy dependence) but not 
necessarily local. Within this schematic model, the diagrams of Fig. [j] are now described by those in Fig. g. 

At t = — oo the state consists of a single incoming pion. In this model the number of pions is conserved by all pieces 
of the Hamiltonian. On the other hand, Hf preserves the number of ph's but F' and F act creating and deleting one 
ph, respectively. A generic state of the system will thus contain exactly one pion plus a number k of ph particles, 
with < k < A where A is the mass number of the nucleus 



hA)=]TMph) fc }:=5>) 



(5.4) 



k=0 



k=0 



| 7t) = \k = 0) corresponds to the elastic channel, and |fc) is obtained after k inelastic steps. Note that in this model 
the ph does not have self-energy graphs (unless absorption is included by a vertex 7rph as in Appendix g). The 
conservation of the number of pions and the free propagation of the ph's are both consequences of the fact that the 
model does not implement crossing symmetry for the pions (there are no anti-pions in this model). In this sense it 



has some resemblance with the Lee model 



52 



where the number of possible graphs is severely limited due to 
the existence of very restrictive conservation laws on the number of particles. In our case the number of graphs is, 
however, considerably larger since the number of ph's is not restricted: in the one-pion subspace the more general 
graph consists of the continuous pion line with zero or more outgoing ph lines stemming from it plus zero or more 
ph internal lines with the ph emitted and reabsorbed by the pion in any order. Since, within the model, a ph cannot 
couple to a pion-anti-pion pair there are no ph self-energy graphs. This model is devised to describe, in a simplified 
manner, the problem of pion-nucleus reactions at energies around resonance or below where pion production is below 
threshold or barely so. When we introduce absorption in Appendix ^ the pion will couple to a single ph (modeling 
the absorption of a virtual pion by a nucleon). In this case pions and ph can be transmuted into each other and the 
number of pions needs not be conserved. (Nevertheless in this region of energies, states with two or more pions can 
only be virtual and we shall simplify further the exposition by not including them in the formulas.) 
Our starting point is the Schrodinger equation 



H\r/>,t) 



(5.5) 
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Using that the states \k, t) are linearly independent, we can write 

ih^-\k,t) =Hf\k,t) + F^\k-l,t) + F\k + l ) t), k = 0, 
at 



(5.6) 



(Here and in what follows we use the convention that quantities with indices k out of the physical range < k < A 
vanish identically. In the present case the terms with \k = —1) and \k = A+l) in the equations for k = and k = A, 
respectively, are absent.) Although correct, these equations show an unwanted symmetry under time reversal: as the 
pion goes scattering through the nucleus, the number of ph produced will increase, then we would prefer that states 
with higher number of ph's had lower ones as a sou rce a nd not conversely. Let us see how to achieve this in the simple 
case of A = 2. Using an energy representation for (pit): 



E\0,E) 
E\l,E) 
E\2,E) 



H f \Q,E)+F\l,E) 

H f \l,E)+F^\0,E)+F\2,E) 

H f \2,E)+F^\l,E). 



(5.7) 



From the last equation |2, E) = \ E — Hf + irjj F^\l, E). The time reversal symmetry is broken by choosing +ir]. 
Substituting in the second equation and using the same method in the resulting equation for |1, E), we end up with 



E\0,E) = Hf\0,E)+F\E-H f -F(E-Hf + iri) F f 

E\1,E) = Hf\l, E) + F (e — Hf + irj\ * F*\l, E) + F^\0, E) 
E\2,E) = H f \2,E)+F^\l,E). 



-l „ 
iti) F^\0,E) 



(5.8) 



The new equations have the desired form; the equation of motion of involves only the states \k) and \k — 1). In 
general the equations are 



with 



E\k,E) = H k (E)\k,E)+F^\k- 1,E) 



H k (E) = H f +FG k+1 (E)F^ 



G k (E) = [E-H k (E)+ir, 



(5.9) 

(5.10a) 
(5.10b) 



Comparing (5.9) with (4.2), we can see that the operator H k (E) plays the role of an optical Hamiltonian for the state 
7r(ph) ). The operators G k (E) are the corresponding propagators and in this model they do not contain intermediates 
states with less than k ph particles. From our previous convention H k=J ^(E) = Hf and G k —A(E) is just the free 
propagator. On the other hand, G k= p(E ) is the full pion propagator and H k= o{E) is th e fu ll pion optical potential. 
For this elastic channel the last term in (5.9) is abse nt and this equation coincides with (4.2). 

Upon Fourier transform in time-energy Eq. ( |5.9|) provides the differential-like evolution equation. It can also be 
written in integral form as 



\k,E) = G k (E)Ft\k-l,E) (fc>0). 

Both forms of the evolution equation will be used subsequently. 
To proceed, we define the density matrices 

p k (E u E 2 ) := \k,E 1 )(k,E 2 \ 



which satisfy (using ( 5.1l| )) 



Pk 



{E X ,E 2 ) = G k {E x )F^p k ^ (Eh., E2)FG\{E 2 ) {k > 0) 



(5.11) 



(5.12) 



(5.13) 



In Wigner's representation, this equation gives p k (t,E) in terms of p k - X (t, E). In order to obtain the differential-like 
equation for p k (t, E) we first rewrite Eq. ( 5.13 ) as 



Eip k (Ei, E 2 ) = H k (Ei)p k (Ei, E 2 ) + F^ p k -i(E x , E 2 )FGl(E 2 ) 



(5.14) 
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and with the same notation and method used to obtain fl4.12p and ( 4.13| ), we have 



up k (u>,E) = H k (E + \u)p k (u,E)-p k (u,E^Hl{E-\w) 

-G k (E + iuj^pk.xF + pk-\F& k (E - \u) , (5.15) 



or in time representation 



ihd t p k (t,E) = e^ P)ai E H) H k (E)p k (t,E) -e- 1 ^ 9 ^ p k {t,E)Hl{E) 



_ e iinay>d^ G k {E)F^ p k _ 1 (t,E)F + e~? lhd > [ 9 * F^p k -i(t, E)F& k (E) . (5.16) 

This can be regarde d as a n extension of Eq. ( |4.13| ) to account for the creation of the new \k) states out of \k — 1) 
when k > 0. Using ( |5.10a ) it is readily verified that this equation preserves unitarity, namely, 



dtj^tr (p k (t, E)) = 0. (5.17) 

k=0 

To summarize this Section: we have been able to write an evolution equation for the quantities p k (t, E), namely Eq. 



(5.16), which in full Wigner's form provides us with the relationship between simulation-like quantities, p k (x,p; i, E), 
on one hand, and microscopic- like quantities, the Green's functions G k (E) and optical Hamiltonians H k (E), on the 
other, without semiclassical approximations involved. However, this equation is not fully satisfactory because too 
much information is contained in p k , namely, all nuclear degrees of freedom as well as far from classical (highly 
virtual) pionic degrees of freedom. We deal with such problems in next Sections. 

VI. REMOVAL OF PH DEGREES OF FREEDOM 

The quantities p k (t,E) and their evolution equations contain information both on "pions" and on " ph " states. If 
we are only interested in the pionic reactions it is convenient to simplify the problem by just working with the pionic 
degrees of freedom. This is the approach in p4[ , where the simulation only traces the path of the pion inside of the 
nucleus. This requires to eliminate the ph degrees of freedom. 

The first idea is to defi ne a new density matrix out of p k for the pion only by taking trace over the ph part. If we 
attempt to do so in Eq. ( 5.1 3| ) , or in the other equations, we find that this new operator tr p h (p k ) does not obey an 



autonomous set of equations; the knowledge of tr p h (p k -i) docs not provide us with tr p h (p k ) because this information 
is only partial. On the other hand, tr p h does not remove the ph energies which are included in E, and in addition, 
the time t is a common time for the pion and the ph's and that may not be the most appropriate choice. It seems 
thus necessary to disentangle the different energy and time dependences in p k in order to find a density matrix truly 
depending only on pionic energy and time, as well as and p^. Likely, the problem is not trivial, and in fact I have 
only partially succeeded in solving it: there is a solution if the propagators G k {E) only contain direct graphs, as those 
in Fig. ||a, and no crossed terms, Fig. |^&. In this case the ph particles are distinguishable and 1, 2, . . . , k labels the 
order in which the ph's have been produced. 

When the particles are distinguishable, the total wave function at time t is just the product of individual wave 
functions taken at the common time t. In energy this corresponds to a convolution. Mathematically we have 

ihG k (E) = [ihGl k) o ihG^ o ifiG^ o • • • o ihG^} (E) + (crossed terms) (6.1) 

where the explicit term is that associated to direct graphs, and we have introduced the following notation: the symbol 
"o" stands for convolution over the energy dependence 

[A o B](E) = J ^ ^27rhS(E -E A - E B )A{E A )B{E B ) , (6.2) 

the operator G^ (E) is the free propagator of the £-th ph (as noted previously, the ph particles do not have self-energy 
graphs in this model) 

G$(E)=(E-H$+ir)Y 1 , (6.3) 
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FIG. 3: a) Typical graph contributing to the directed part of Gk', b) crossed graphs in Gk (for k = 1). 



{Hph i s the free Hamiltonian of the £-th ph, included in j?Q Ph - ) ). Finally, the operator G^(E) is the pion propagator 
in presence of k ph's, but including direct self-energy graphs only. Thus it is given recursively by (cf. Fig. |^a) 



G^{E) = (e-H^{E) 



IT) 



H£\E) = -inP[G^ odi k h +1) ](E)F^ 



(6.4a) 
(6.4b) 



where (E) is the optical Hamiltonian of the pion in presence of k ph's and the operators F, F^ act only on the 
k + 1-th ph. Note that G^] and G„ act in the £-th ph and pionic Hilbert spaces respectively and they commute. 

The purpose of selecting the direct graphs was to be able to disentangle the energy-time and position-momentum 
degrees of freedom carried by each individual particle in the state p k (t,E). This is achieved as follows. As shown 
in detail in Appendix |B], if crossed terms are dropped in Gk, a new (more detailed) state \E w ,Ei, . . .,E k ) k can be 
defined, depending on ir and ph degrees of freedom (a; and p) and also on their energies, which is related to the state 
\k,E) in Eq. (|5.9|) by means of 



(6.5) 



\k,E)= [ ^^±...^2irn6(E-E„-E 1 E k ) \E«, E u . . . , E k ) h . 

J 2nn 2ttTl 2-kJi 

The new state satisfies the integral-like equation, which is analogous to ( |5.11 ), 



\K,Ei, E k )k = ih&V t \E*)ih&$(E k )±:P*\E* + E k ,E!,E 2 , E k -i) k -i (k > 0) . 



(6.6) 



Here F f acts on the pion in | ) k -i and creates the fc-th ph, thus building the subset of direct graphs only (and so 
Ek,Ei, ... ,Ek)k is not symmetric under exchange of ph's). Actually this equation is used recursively in Appendix 
to define the states \E n , Ei, . . . , E k ) k starting from the k = state. 
Next we define the associated density matrix 



Pk(^TT, E n ,u)i,Ei, . . . , cj k , Ek) 



\E 7T 



5^tt, E\ + . . . , E k + \uk) k (En — Ei — \tjii, ... ,E k — ^WfeL 



(6.7) 



from which p k {to, E) can easily be recovered by making use of Eq. (6.5). This new density matrix satisfies the following 
recurrence equation (which is a translation of (|6.6|)) 



p k (^,K, oji, Ei, . . . , Wfc , E k ) = h 2 Gi k) (E„ + i^) (E k + |w fc ) 

xp fe _iK +uj k ,E„ + E k ,uji, Ei, . . . ,Lu k ^i, E k -i) F G ( £HE k - ±c^)Gl fe)t (K - ^) . (6.8) 



15 



This equation represents a definite improvement over Eq. (5.13) because now the trace can be taken over ph de- 
grees of freedom and a closed set of equations is obtained: on the right-hand side the trace factorizes in the form 
{1, . . . , k — 1} {k} and as a consequence tr p hPk is given in terms of tr p hPk-i- 

In order to obtain a pionic density matrix, let us call it p^ k '(u>n, En', E), the ph energies Eg can easily be integrated 
out. Also, the "times" uig should be fixed. After Fourier transforming, each cog becomes the time tg at which the £-th 
ph is detected. In order to preserve unitarity at any "pionic" time t, we should have t smaller that any ph time, so 
we are lead to choose all the ph times tg as +oo in the definition of the purely pionic density matrix, that is, 



p^(t,E^E) 



lim 

{t £ }^ + oo 



dE x dE 2 dE k 



2irh 2ttH 2irh 
xtr p hPfe(i, Ek] t\,Ei, . . . ,tk,Ek 



2nh S(E -Ei 



E k ) 



(6.9) 



(Where t refers to the pionic time associated to the pionic frequency w,. Also we remark that we are taking the limit 
of large ph times, and not integrating over those times.) This is the density matrix in pionic space describing the 
pions which have scattered k times (producing k ph's). In particular, 



p (k ~ Q) (i, En, E) = 2irh 6{En - E) p k=0 {t, E) 



(6.10) 



is the density matrix for the clastic channel. Clearly, in p^ k \t, E v ; E), E corresponds to the total energy, which is 
conserved. On the other hand E v is the energy carried by the pion after k collisions and both quantities coincide 
in the elastic channel. We further discuss on the coexistence of these two energies in p^ (t, E v ; E) over the end of 
Section [VD. 

As shown in detail i n Ap pendix |c], if the definition in ( |6.9| ) is used in Eq. ( |6.8[ ), the following recurrence is found 
for p^ (analogous to ( 5.13 )): 



J ^^{ 2 ^^-E n -H ph )F^^- 1 Hu J ^E' 7r ;E)F} 



xG^iEn-iun), (fc>0) 



(6.11) 



Moreover, there it is also shown that unitarity is preserved, as a direct mathematical consequence of taking all i p h as 
+oo, that is, 



A 

E tr 

,fe=0 



dEn .( fe 



2nh 



= 0. 



(6.12) 



To obtain the Schrodinger-like equ ation associated to the integral- like equation in ( 6.11 ), we follow the same 
procedure as that used to obtain Eq. (5.15). This gives 

^K,£,;£) = H (k) {E^ + \^)p {k) {^E^E)-p^{^,En-,E)Hn k) \En - §ug 

/dE r 
-^27rhS(E^-En-H ph ) 
Znfi 



X G^{En + \^)F^ p^){Un,E'^E)F 

-Ft p(k-i) F4; F) F Gi fe)t {En - \Un) 
and then, after Fourier transforming in u>n, 

d t pW(t,En;E) = (d t pW) (+ \t,En;E) + (d t pW) ( ~\t,En-,E) 



(6.13) 



(6.14) 



with 



^hd t p^) i+ \t,En,E) = e^ s ^Hn k \En)pW[t,En-,E) 

-e-i ihd ^ d ^p( k )(t,En,E)H^(En) 



(6.15) 
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FIG. 4: Quasielastic process in a finite system. 



and 



(ih d-, 



hp 



(*) 



(-) 



(t, E n ; E) 



-tr 



2?r?i 



(6.16) 



Equations ( 5.1l| ) and (6.13) (or equivalently, ( |6 .14 )) are the relevant result of this Section. In dtp^ we have dis- 
tinguished two contributions. The first one, {dtffi k ')( + \ is related to p^ itself and describes a pion (of class k) 
propagating with optical Hamiltonian H^ k \ (cf. Eq. ( 4.13 )): it contains both, the "elastic" propagation of the pion 
inside the nucleus, and the quasielastic steps in which the pion of class k becomes of class k + 1. In this sense it can 
be called the annihilation part of dtp^ ■ The second contribution, (dtp^)(~\ accounts for the quasielastic steps of 
the form k — 1 — > k. From the point of view of it is the creation part of dtp^ ■ The Dirac delta in this term 
indicates that the ph's are on-shell. Mathe mati cally, this is a direct consequence of having chosen all the ph times as 
+oo in the definition of the pionic density, (6.9). 

Before studying Eq. ( |6.14 ) in more detail, let us see how it works in a simple case: 



F= I d 3 xg(x)$ pll (x)4>l(x)^(x) (6.17) 

and Hq , .Hph functions of the momentum only. Let us first put the creation part of dtp^ in Wigner's form. In a 
first step 



/HE 1 
— £ 2nh 8{E'„ -E w - H ph ) p^ (t, E'^E)F 
Znn 



where 



Vhkf (2^)3 ^ ? ( x *> p * + P' + * *> E * + H ph ( P ');E) 



g 2 (x,q) := J d 3 ye-*yV h g*(x + ±y)g(x-±y) 



(6.18) 



(6.19) 



and q is the momentum transferred to the nucleus (Fig. |j). 

Next, we shall retain only the leading order in Ti. In this way several simplifications take place in (d t pY^'- i) 

the exponentials exp(±^ih~dj. p ^ d^) are unity at leading order in h; ii) the pion propagator becomes free, G^iE) — > 



(E — Hq + irf) 1 , because e ach loop in its self-energy gives a fi, factor (Eq. ( [4.16 )); iii) in the classical limit the 
operators commute (cf. (Eq. (3.S)); and, iv) using the identity (3.5) 



g 2 (x,q) -^ (2nhy S(q)g 2 (x) 



(6.20) 
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After some algebra, the creation part can be written as 

(-) 



with 



Q(x, p, E n ; x',p', K) p^ Or V; t, K;E) 
(2irh) 2-kTi 



O(h) (6.21) 



Q(x,p,E^x',p',E' v ) = -^g 2 (x) 5(x - x')2n?i5(E n - H^'(p))2TTnS(E' n - E n - H ph (p' - p)) . 

n 



(6.22) 



We can see that the quasielastic probability in this approximation is positive, local (x = x') and momentum is 
conserved (q = 0). The ph are on-shell (E^ — E n = iJ p h(p' — p)) and moreover the outgoing pions are also created 
on-shell. As noted in the Introduction the latter fact is a typical consequence of taking ft — > 0, due to the uncertainty 
principle. Mathematically, the factor 2irTi5(E 7V — i/^*' "*(£>)) follows from the imaginary part of the free pion propagator 

G^iE) — G\{E) in ( |6.16 ). Summarizing, the probability of quasielastic Q obtained at leading order in ft coincides 
with the result that would follow from carrying out a standard nuclear matter calculation (g = g c constant) in lowest 
order in perturbation theory, plus a local density prescription (<? c — > g(x) at the end of the nuclear matter calculation) 
pjj . This latter fact is quite remarkable since usually the local density approximation is put in by hand, and here it 
follows naturally as a semiclassical approximation. 

The annihilation part of dtp^ comes from computing the pion optical potential up to one loop and keeping the 
free and the imaginary parts (lowest order in ft in Eq. (4.17)). This gives 



(d t p {k) ) i+) (x,p;t,E^;E) = {i#°(p), P (k) (x,p;t, E„; E)} - R(x,p; E„) p (k) (x,p;t, E„- E) + 0(ft) (6.23) 



where 



R(x,p;E 7r ) = -^g 2 (x) I tohSfa - H ph (p -p') - H^ip')) 

n J (27rnj 



d 3 p' 



(6.24) 



The first term of the right-hand side of Eq. ( |6.23| ) describes a classical free propagation of the pion. The second term 
indicates a reaction probability rate given by R{x,p, E„). Unitarity is verified since 



R{x,p;En) = 



d 3 x'd 3 p' dE'„ ~ , , 
(2irh) 2nh 



(6.25) 



Of course, the right-hand side of (6.21) should be set to for the elastic channel, k = 0, and similarly, the reaction 
term is not present in ( 6.23 ) when no ph remains in the nucleus, k = A. Otherwise, Q and R are independent of k 
and E in t his ap prox imation. Also note that bec ause the pions are on-shell the E n dependence can be dropped and 
Eqs. ( 6.21 ) and ( |3.23| ) have the form of Eq. ( |l.3| ), except that the classical propagation has been extracted from the 
quasielastic rate Q. 



Unfortunately, beyond leading order in ft, Eq. ( 6.14 ) displays two undesirable features: first, the dependence of 
p^ on En, which implies that virtual as well as real pions coexist in p( k \ We deal with this problem in next Section. 
And second, the non locality in time: in general dtp^ k ' will depend on the previous history of and p( fe -i) . The 
instantaneous equations can be obtained by a method similar to that used for Eq. (4.17): dtp^ is given as a function 
of d?pW an d a™'p( fc - 1 ), n,n' — 0,1,..., each time derivative carrying a factor ft. Then, hi gher order derivatives 
on the right-hand side can be eliminated in terms of lower ones. For instance, writing Eq. (6.14) with an obvious 
schematic notation, we have for k = 1, 2 



d t p 



(2) 



(1) 



(6.26a) 



Then, to first order in ft, 



ftpw = y: a- n>i +) ap p w + y: ft" Nt ] s t y o) 



n=0 



{hN ( - ] N' { - ] + 0{Ti 2 )} p 



(0) 



(6.26b) 



(6.27) 
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We observe that the equations can be written in an instantaneous form at the price of expressing dtp^ , not only in 
terms of p^ and p^^ 1 ^, but using all densities p( k ' with kl < k. We shall see that a similar phenomenon occurs 
when the E v dependence is removed form p^ k \ 



VII. INTEGRATION OF THE VIRTUAL PIONS 



In the previous Section we have succeeded in writing an equation, Eq. ( 6.14 ), for the temporal evolution of the 
density matrix of pions which have scattered k times, with only pionic degrees of freedom. However, it is not a fully 
satisfactory one due to its dependence in the pion energy, E v . Consider, for instance, the evolution of a pion after 
its creation (i.e. after the scattering): two frequencies are relevant, first that characteristic of the source, represented 
by E n , and second that of the free evolution, H^\p^). Due to interference, after traveling some wavelengths only 

on-shell pions, those with E v — Hq(p v ), will survive, unless another scattering takes place. Such a collapse of the 
wave function can also be regarded as a classical limit: if h~ were very small any time interval would be large compared 
with a few periods of the order of H/\E V — H^(jp v )\ provided that the particle is off-shell. Then we expect the 
following relation to hold: 

pW(t,E v ;E) -^2wfiS(E 7r -H^) (7.1) 

ft.— j-G 



up to a fa ctor. This relation follows from the energy-shell constraint Eq. (4.21) and is checked in next Section (cf. 



Eq. (|8.13[)). 

Initially (in the incoming pion beam) the pion is on-shell. After several quasielastic steps the pion leaves the nucleus 
and is again on-shell since only on-shell pions can travel long distances to the detector. Between two successive 
collisions the pion can be (nearly) on-shell (i.e., a real pion) or off-shell (a virtual pion). Consider a typical path of 
the pion through the nucleus, for instance RVVRVR (time increasing from right to left) where R stands for real and 
V for virtual pion, in this case with five quasielastic collisions. Each collision is of the type one-body mechanism, 
producing one ph. The virtual pions will travel a short distance and the Monte Carlo simulation only needs to trace 
the real ones. Within a classical limit, the real pions would follow a classical trajectory, whereas the path followed 
by the virtual pions shrinks. The situation is thus better described as a path RRR, that is, involving real pions only 
and two collisions. The first collision ejects two ph, a two-body quasielastic mechanism, whereas the second collision 
involves a three-body quasielastic mechanism (see Fig. ||). 

Advancing results of the next Se ction , we can see the different properties of real and virtual pions as follows. Consider 
the classical limit directly in Eq. ( |6.11 ). By neglecting the ±o> w /2 (which carries an K upon Fourier transform to time 
representation) in G v and going to time representation it turns out that p^ (t, E v ; E) is related to p( k ~^ (t, E' v ; E), 
i.e., the number of pions of class k is proportional to the number of pions of class k — 1 at the same time. This is only 
possible if they live during a very short period of time, and this is correct for virtual pions. Of course, for on-shell 
pions this would be an absurd consequence which is avoided because for them taking ±u> n /2 — > is not justified even 
in the classical limit: the quadratic G^G^ diverge nce (a t the on-shell point) is too strong. (This is also the reaso n 
for needing differential-like evolution equations, as ( |6.15| ) which have softer divergences than the integral-like (6.11).) 



The virtual pions are not only irrelevant to final cross sections but their short living times can only be achieved by 
interference which is very hard to reproduce by a simulation. Also, their existence implies that p(t, E„; E) should be 
strongly non-positive definite. 

Our goal is thus to integrate out the virtual pions and end up with an equation including only real pions. Because 
the real pions are on-shcll their energy dependence is fixed and its corresponding density, to be denoted fk(t, E), does 
not depend on E v . In addition, in order to preserve unitarity at each time, and not only in the t — > +oo limit, we are 
led to the following definition for the density of on-shell pions, 

/*(*, E) := J (t,E^;E) , k = 0,1,2,..., (7.2) 

as the Monte Carlo density to work with. The suitability of this definition is further discussed over the end of this 



Section. By construction it satisfies a conservation equation like Eq. (3.12), namely, 

A 



^2tvf k (t,E) 



.k=0 



= 0. (7.3) 



It is clear that the evolution equation for p^ h \t, E n ; E) completely determines the evolution of fk(t,E), however, 
whereas the equation for p^ involves pions of type k and fc — 1 (one-body mechanism) that for /& will depend on 
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pions of all classes k' < k (A^-body quasielastic). In addition, initially p^ is determined by f (cf. ( 6.1 0| ) ) . So there is 
an autonomous equation for the densities The question arises how to obtain such an equation. Simple integration 
over En in Eqs. (3.11) or fl6.13| ) does not work. Essentially the problem is how to invert Eq. ([7^). To do so I make 
an Ansatz for the E v distribution of the real pions in p^ k '(t, E n ; E): 



pF> (t, E„;E) = p (k) (t, E 7r ;E) + (t, E^ E) 



(7.4) 



where pr.v are the distributions of real and virtual pions in p. The explicit form of these distributions is restricted 
by imposing the following conditions: 

1) Pft\t, E„;E) should be constructed out of fk{t, E), in order to be able to invert Eq. (7.2). 

2) Consistency with Eq. ([7.2[) requires that 



fk(t,E) = 
= 



dE, 



7T Jk) 



2nh 



p%>(t,E w ;E) 



dE * -(fe) n P .m 



(7.5a) 
(7.5b) 



because fk is already the distribution of real pions without the redundant E n dependence. 

3) p^ , py 1 should be Hermitian. 

4) In the classical limit, py should collapse to zero, due to quantal interference: 

P P(t,E^E) — ♦ 

pf(t,E v ;E) — > 2TrhS(E 7r -H i * ) )f k {t,E). 



(7.6a) 
(7.6b) 



The (S-function in Eq. (7.6b) follows from (7.1) and the f actor / is obtained by normalization. In the classical 
limit, operators commute so there is no conflict between (7.6b) and the point 3) above. 



5) For the elastic channel we impose the constraint 



pt 0) = o. 



(7.7) 



This defines the elastic channel pions as real. Virtual pions only appear as a consequence of (hard) collisions, 
whereas the mean field effects under the elastic evolution are regarded as soft. Technically this choice is needed 
to be able to close the equations below. 

Let us see how such an Ansatz solves the problem. To alleviate the notation, let us write Eqs. ( 6.11 ), ( |6. 13 ) in the 
form 



M (k,k-1) p(k-l) 



or even more compactly 



p = Mp 
dtp = Np 



(7.8a) 
(7.8b) 



(7.9a) 
(7.9b) 



where p is a (column) vector (p'- *' , p^\ . . . , p^) , and M,N are m atric es with respect to the index k and super- 
operators in Hilbert space (map operators onto operators). Also, Eq. (|7.2|) can be written as 



f = bp, 



dEn 
2^h ' 



In addition, p R will be constructed out of fk by some linear procedure to be specified later: 

PR = af = abp , 

pv = (1 - ab)p := P v p. 



(7.10) 



(7.11a) 
(7.11b) 
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>0) 



i is the projector onto py\ Equations (7.5) then read, 



ba = 1 . 



bP v = . 



(7.12) 



Note tha t b removes dependence on E v while a creates dependence on E^, and M increases the index k by one unit. 
Eq. (7^4) then implies 



PR + Pv =af+ P v p = af + P v Mp . 



(7.13) 



In the last step we have used the integral-like evolution equation for expressing the short-living virtual state in terms 
of its source. (In the elastic channel this is consistent due to fl7.7| ).) This process can be iterated until the source is a 
real pion and in this way the virtual pions are integrated out. In practice this procedure can be carried out as follows. 
Using previous equations and the fact that dt and b commute with each other, we can write the following chain of 
relations: 

d t f = bd tP = bNp 

= bN (af + P v Mp) 

= bN a f + bN P v M (a f + P v Mp) 

= bNaf + bNP v Maf + bNP v MP v M(af + P v Mp) 

= ••• . (7.14) 

At each step, py is carried to smaller values of k since it appears with a new power of M. This downward recurrence 
ends due to (7.7) (the elastic channel pions are real and no inconsistency arises). In this way all virtual pions can be 
eliminated and formally 



d t f = bN(l-P v M)- 1 af. 



(7.15) 



This equation explicitly shows that /, as a vector, satisfies a closed equation. Let us expand Eq. (7.15), with obvious 
notation, 



dtfk 



{&JV(M) a (*0} f k 

+ | 57V (fe,fc-i) a (fe-i) +feiV (M)pW M (M-i) a (/c-i)| j k i 

+bN {k,k) p p M {k,k-l)p£-l) M{k -X,k-Z) a {k-2)\ f k 



dtf 



R {k) h + Q {k ' k - l) h-x 

Rf + Qf 



Q 



(k,k-2) 



fk-2 + 



(7.16a) 
(7.16b) 



where R tak es c are of the propagation of the pion once it is produced and Q describes the quasielastic steps. (Note 
that in Eq. (1.3) the elastic channel was included in Q while here it is given by i?( fc =o) ) 

Our actual procedure has been to remove any virtual pion by explicitly writing it in terms of the nearest real pion 
acting as a source for it. As a consequence only real pions do explicitly appear, and the hidden virtual intermediate 
states translate into effective iV-body quasielastic probabilities, Fig. [|. 

I n a cla ssical limit, each factor Py will increase the order in fi, due to Eq. ( 7.6a| ), and only few terms will be needed 
in ( 7.16a ). For instance, at leading order in h the one-body mechanism, Q^ k ~ 1 \ is dominant and the effective 
two-body quasielastic, Q( k ' k ~ 2 \ is the first quantum correction. 

R^ plays the role of an effective Hamiltonian for //.. Its actual form will depend on the concrete form of a^- k \ 
which still has to be chosen, but using general properties of we still can say something about R^: 

i?W/ fc K,S) = (bN^a^f k ) K,£) 

= ^ / ^(H£\E r + \u,Jj%\u,„E*^-pf 



1 

ift 



, fk {^TT , E) 



J S (^ fc(K + h u *)$lHu*,^E)-$l\w„E„;E)v2(E„ - ^)) (7.17) 



1 

ih I 2Trh 
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R 



FIG. 5: Shrinking of virtual intermediate states (V) to produce an effective iV-body quasielastic step between real pions (R). 



where Vk — — H^ 1 is the optical potential. Further, for fi — > 0, using Eq. (7.6b) 



iTi 



H^iK + \^) f k (LJ n ,E) - / fc K, E) H£>\K 



(7.18) 



where E v = Hq is the on-shell energy. 



In the special case k = 0, Eq. (7.7) completely determines a^ -0 - 1 with the help of Eq. (6.10) 

a ( fc =°) =2nh S(E-E^) 



(7.19) 



whic h obv iously satisfies all other requirements 1) - 4), and the equation dtfk=Q — R^ k ~ ^ fk=o is nothing else than 
Eq. ( pL13l ). 

In the general case, the super-operator a^ k ' can be chosen fairly arbitrarily. A natural choice is 

(7.20) 



(7.21) 



where U^if) is the evolution operator, and unitarity, i.e. U v {t = 0) = 1. If ?i ^ 0, G v and / commute and G^ 
approaches the free propagator, 



(a^jk) (t, E^E) = in&W(E«) / fc (t, E) - f k (t, E) in&™(E K ) , k ? . 
Clearly, this choice satisfies the requirements 1) and 3), above. The point 2) follows from 

ihG 7r {E 7T ) = I ' dte lE ^ /h U^{t)9{t) 



mG^{E v ) - ih&W(E„) — > 2irhS(E v - H^>) 



(7.22) 



which is the point 4), above. Other choices are possible, for instance taking only ImG,r, or taking G* cc . Another 



possibility would be to use E v ± |w T instead of E v in G^ in (7.20) (in frequency representation), however this choice 
turns out to be inappropriate for computing iV-body absorption processes (cf. next Section). As compared to that, 
our choice of is the instantaneous version since it relates to fk at the same time. It is noteworthy that the 
definitio n of has some resemblance with the formulas invoked in the standard cascade method jj], O. 

Eqs. ( 7.16| ) are, perhaps, the main result of the paper. They describe the evolution of the purely pionic density 
matrix containing only real pions. They make explicit the iV-body quasielastic rates Q^ k ' k > seen by the pion as it 
propagates through the nucleus. These rates are the input to be used in the cascade method. In Appendix ^ these 
formulas arc trivially extended to include absorption. They are analyzed in next Section. 

It is noteworthy that the concrete choice of cann ot a ffect the evolution of fk, if computed to all orders, that is, 
including all iV-body mechanisms, since its definition (7.2) is independent of a^ k \ (The right hand side of (7.15) is 



actually independent of a.) Obviously, different choices introduce different organizations of the series. This ambiguity 
is related to that in the separation into real and virtual pions. In this regard, it would be very interesting to choose a^ k ' 
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so that not only bP^ 



(fc) 



0, but also 6iV( fe ' fe )p| fe) 



0. This essentially means that the real and virtual components of 



the density are separately preserved under elastic evolution, i.e. in the absence of collisions. In this case Eqs. (7.16a) 
would simplify considerably. However, it is not clear how to impose this property or even whether it is consistent 
with the other requirements set o n . 

We also note that the formulas ( 7.16 ) involve no approximations, except that of not including crossed graphs, that 
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is, neglect ing Bose symmetry of the ph's, as explained at the beginning of Section 

In Eq. (|6.9|) we introduced the density matrix p^ k \t, E^; E) for class k pions, depending on the time and e ne rgy o f 
the pion a nd al so on the conserved total energy E. It can be noted that in most f ormu las, including ( |6.11 ), (3.13), 
(7.2) and ( 7.20 ), E appears only as a parameter in the densities. An exception is ( 7.1S ), which refers to the elastic 
channel, and so a separated discussion is needed for that case. 

For pions of class k > 1, E is just a parameter. We could as well introduce a new density p^ (t, E^), without 
E dependence, by integrating p( k > (t, E^; E) over E, and rederive all equations for it. If this is done, the density of 

corresponds to the standard 



real pions defined in (7.2) becomes fk(t), which, recalling the "energy average" in 
equal-time pionic density for pions of class k. In this light the integration of virtual pions is related to the closing of 
a set of equations for the equal-time densities. After the pion has scattered once or more, the initial energy E is no 
longer relevant and the corresponding A-body reaction probabilities do not depend on E, a fact explicitly verified in 
next Section. 

On the other hand, for pions in the elastic channel E is quite relevant and coincides with E„ (cf. (6.1C )). Integrating 



over E„ to obtain fk(t,E) in (7.2) merely removes this redundant energy dependence. It should be noted that, being 
conserved, E is a known datum in the collision experiment, thus it does not seems advisable to remove this information 
by going further and integrate over E to work with the equal-time elastic channel density matrix. Indeed, as verified in 
next Section, the A-body reaction probabilities Q( N ' QS> and A^- N ' QS> display an explicit dependence on the pion energy. 

As noted in Section [V, when a space-time Wigner form is used (as opposed to an equal-time formulation) it is 
customary to use a transport and a energy-shell constraint equation. In this work we have used the differential- 
like transport equation, ( |6.13 ), and the integral-like equation (6.11) as an equivalent set of equations. Actually the 
integral-like equation contains both the transport and the constraint equations (except for the elastic channel) and 
not surp rising ly this equation is responsible for putting the particles on-shell as fi goes to zero in the present approach 
(cf. Eq (§13)). 



VIII. N-BODY EFFECTIVE QUASIELASTIC AND ABSORPTION MECHANISMS 



In this Section we shall work out the consequences of t he p revious scheme in the simplest ca ses. To simplify we 
shall study only nuclear matter (g(x) = constant = g in ( |6.17 ), K(x,y) — constant = n in (Dl), and p 7r (x,p;t, E) 
independent of x), free propagators, locality in time, and in general lowest orders in h. As was shown in Section VI 
all these assumptions are compatible with fi — > 0. 

Let us study first the one-body quasielastic Q^ 1 ' ); 



Q (i,o) = 6A r(i.o) a (o) + & AT(i.i) P (i) M (i,o) a (o) 



(8.1) 



In leading order of ft, only the first term contributes (besides in nuclear matter and using free Hamiltonians N^ k ' k ^ = 0). 
The result can be obtained more directly from Eqs. fl^2l| ), ( |6~2^ ): 



<5 (1,0) (x, p; x', p'; E) = ^ g 2 S(x - x') 2nh S(E - (p) - H ph (p' - p)) 



•2) 



The one-body absorption rate, A^ 1 ' ) is (the operation P plays a similar role in absorption as N does in quasielastic, 
see Appendix O for details) 



that can be worked out to give 



A (i,o) = Po (0) 



•3) 



A< x '°> {x,p; E) = -2 2-kTi S(E - H ph (p)) . 
n 



(8.4) 



However, in leading order in h, E = Hq (p) what is incompatible with the (5-function in Eq. (3.4) in physical cases: 
the pion and the ph cannot both be on-shell and as a consequence there is no one-body absorption in the classical 
limit. 
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FIG. 6: a) Quasielastic step: the incoming pion and outgoing ph are on-shell. The outgoing pion can be real or virtual. 6) 
Absorption process. 



In order to study more complicated cases, a drastic simplification in the notation is convenient. Consider the 
relation = M^ 1 ' )^ ) /o which is represented in Fig. ||a. Using ( |6.11 ) for M^ 1 ' ) and ( 7.19| ) for a^°\ it can be 
written as 



fP-\u*,E K ;E) = G^(E n + ±c^)tr ph {2nh5(E - E„ - H ph )pi / K, E)P} G^(E^ ~ . 



(8.5) 



In what follows, we will use E to denote the total energy E and fo(t) to denote (the full Wigner's form of) /o(t, E ). 
In addition we introduce the notations (see Fig. ||a) 



E := Ev-H^ip) 
E := E Q -H ph (q)-H^\p) 
G(E) := (E + i?])- 1 . 



(8.6) 



Furthermore, we will focus on the fi and energy-time dependence. Integration over momenta, vertex operators (F, 
K), traces, etc, will be implicit. Schematically, the relationship p^ 1 ' = M( 1,0 W°) fo will look (in time representation, 
where the product of frequency functions becomes a convolution over times) 

(i, E) = J dr h 2 (GG 1 ) (r, E) f (t - t) 2t:% S(E - E) . (8.7) 

The total energy E (called E in previous Sections) is not displayed. The operator (G&) (r, E) comes from G^E^ + 
^uj w )Gl(E v - |w w ), namely, 



H 2 (G&)(t, E) = J ^e-^ T ' h h 2 G{E + \uj)G\E-\uj). 



(8.8) 



This is the kernel of the super-operator M and controls how (real and virtual) pions of class k are produced out of 
pions of class k — 1 and their subsequent propagation. In terms of a standard diagrammatic calculation (for instance 
in nuclear matter) this would correspond to compute the graphs associated to the transition amplitude between the 
initial and final states and then to square it to get the transition probability, or equivalently, to compute the self- 
energy graphs of the initial state and apply Cutkosky rules to associate each cut of the graph to a contribution to the 
transition probability to a concrete fina l state. 
The exact integration over lu in 



our purposes. To get such an expansion we start with the identity 

h 2 d T (GG 1 ) (r, E) =ih(G~ &) (r, E) 



is easily performed, but an expansion in powers of h is more convenient for 

(8.9) 
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The right-hand side is the kernel of ./V and is a softer distribution than (GG'Y Some algebra yields 



i%(G — &)(r,E) 



-2klm 



5{t + \ihd E ) G{E) 



2irhS(T) 6(E) + h 2 S'(T)P^r + Oth 3 ) 

E A 



(8.10) 



where 



P Ek : =^(G(i?))' 



l) k d*- 1 1 



.11) 



(k - 1)! dE^ 1 E 
is a renormalized principal value [[35) . Integration over r then yields 

h 2 (G&)(T,E) = -2Mm\e(T + ±ihd E )G(E) 

= 2Trh9{r)6(E) + Tl 2 8{t)P^ + 0(h 3 ) . (8.12) 

E l 

Comparing ( B.10 ) and ( |8.12| ) we can see that ih(G — G^), and so the super-operator N, is instantaneous at leading 
order in h whereas h 2 (GG^, or M, is not. The non-instantaneous piece in M would have been missed if we have 
taken a formal classical limit neglecting the terms in fi. 2 (GGt). This is incorrect because the divergence at the 
on-shell pole is not integrable. Correspondingly, as we will see subsequently, the non-instantaneous piece contributes 
only to real pions and not to virtual ones. Substituting in Eq. ( |8.7[ ), 



{t,E) =2irh8{E-E) [ 2nhS(E)d^ 1 + h 2 -^ + 0(Ti 3 ) ) f (t) 



.13) 



where d t 1 = dt, and the symbol P is implicit. We can see now what fx, p^ and py > look like. Recalling 



a 



b 

(fc) 



ih 



dE 
2^h' 

G-& 



(r,E) = 2%hS(E)d(r), fc^O 



(8.14a) 
(8.14b) 



where in aS > only the instantaneous part should be taken due to the missing ±uj 7r /2 in j7i2Q|) , we have (presently aS > 
acts by convolution over r) 



h{t) = \2TrfiS(E)di 1 + n 2 w + 0(h 3 )j f (t) 
p { ^(t 7 E) = 2nh 6(E) (2irh6(E)d^+h 2 -L + 0(h 3 )) ./„(/! 



)\t,E) = 2nh 3 -L(5(E-E)-S(E))f (t) + O(h 4 



.15a) 
.15b) 
;.15c) 



We can see that in p^ (t) the pions are on-shell, E — 0, and also that it contains a piece, that with d t l fo, which 
depends on the whole previous history of the incoming pions, i.e., p^ (t) contains long lasting components indicating 
that the real pion at t could be produced from a quasielastic step which occurred some time ago. On the other hand, 
Py^ does not contain such components; the number of virtual pions depends on the instantaneous number of incoming 
pions indicating that the virtual pions are short-living states. 

Now we can easily compute absorption and quasielastic rates. Q*- 1,0 ' will follow from 



Q (1 ' 0) /o(*) = J —dTth(G-&)(T,E)2TrhS(E-E)f a (t-T) 
= (2irh6(E) + 0(h 2 ))f Q (t) 
which coincides with (8.2). Analogously, A^ 2 ^ and A^- 2 ' ^ 

A^h = PaW/ 1 = J ^2«h5(E-E A )p$(t,E) 
A^fo = P Pv ] = J ^2«h8(E-E A )$>(t,E) 



(8.16) 

;8.17a) 
8.17b) 
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P, E n y 

p\ ,s q> H vh (q) 



p+q, e 



^', H ph (p') 



(t 2) \ / 



\P,E % ' 



fn 



p+q, e 



a) b) 

FIG. 7: a) Two-body quasielastic. 6) Three-body absorption. 



where Ea ■= i? p h(p) — fig (p) (see Fig. 0o). Using the expressions in ( 8 . 1 5| ) and the fact that Ea ^ 0, as discussed 
for A^ x '°\ we get 



A^f (t) = 27rh 3 -^5(E A ~E)f Q (t) + 0(h 5 ) 



(8.18a) 
(8.18b) 



q(!.o) an d ^4(2.°) are then the lowest order quasielastic and absorption mechanisms. Both are positive definite [P-gs 
is positive outside the pole) |6(|. 

To display further features of the present scheme we shall study the coefficients Q^ 2 ' ) and A^- 3 '°\ Again schemati- 
cally 



Q (2A) h(t) = J^dr'^in(G-&)(r',E')2wh6(E-E'+E')p^(t-r',E) 
= (2irh5(E') + 0(h 2 )) f±(t) 

where E' = E' % — H^(p>), E' = H^\p) - H^\p') - H ph (q') (see Fig. 0a). 
Q( 2 >°) is obtained by taking py* instead of p^\ in ( 8.19 ): 



(8.19) 



Q^°~>f (t) = 2^1 I (5{E' +E)- 5{E')) f (t) + 0(h 4 ) . 



.20) 



Note first that Q*- 2 ' 1 -* = Q( 1,a \ at least in leading order. In general Q( k + S > s, > = Q( fe >°) would be most desirable, but 
given the asymmetry in the definitions of a^ fe ^ ^ , and it might not b e true for higher orders in fi. Secondly, note 
that Q^ 2 ' ** is not positive definite. This is a consequence of Eq. ( |7.5b| ), and it was expected because it is the first 
quantum correction to the quasielastic. Indeed 



Q^fi(t) = 2irh5(E')2irh5(E)d^ 1 fo{t) + 0{h 3 ) 



.21) 



which shows that effectively Q( 2 >°) w fiQ^ 2,1 \ It is also interesting that as E — » in (8.2C) and the intermediate pion 
(line p in Fig. 0a) approaches its mass-shell, the two ^-functions cancel: this is the effect of subtracting real pions 
from . Then, even if the P-^i distribution were not already renormalized and finite [j67| the subtraction pP^ — p^ 
will produce a finite result. 
Similarly for absorption: 



A^h = J d ^ns { E>-E>A) 

Ai3 '° )fa = J 2%?^ E '- E ' a) 



(2) 
PVR 



(2) 

Pvv 



..22a) 
.22b) 



26 



with E' A = H ph (p') - H^\p') ^ (see Fig. |&), and the result 

A^f x {t) = 2TTh 3 ^6(E' A ~E , )f 1 (t) + 0(h i ) (8.23a) 
A( 3 -°)/ (i) = 2vh^ ?t(5(E> A -E>-E)-5(E> A ~E>))f (t) + O(h 7 ). (8.23b) 



Once again A^ 3,1 ' = ^4( 2 '°), also A' 3,0 ) is not positive and the 6's cancel as E — > . Furt hermore, in both cases Q' 2 '°) 
and A^ 3,0 ' average to zero (as functions of £"), again a direct consequence of ( 7.5bj ). Although after momentum 
integration Q^ 2,0 ' and A^ 3 ' ^ will not be zero, they will be small. In our formalism this is reflected in the fact that 
both have higher order powers of h. Remember that fi — > should be understood as a physical limit: it is not fi that 
is small, rather the coefficients with higher orders of h are smaller if the classical limit applies. 

Even though Q*- 2,0 - 1 and A^ 3 ' ^ involve the propagation of intermediate (virtual) pions, both rates are instantaneous 
in leading order implying that such states last a short time. A more detailed treatment equally shows that they run 
a short distance, i.e., Q and A are local at leading order. 



IX. SUMMARY AND CONCLUSIONS 



In previous Sections we have reduced a com plex many-body evolution equation to something more similar to the 
master equation of the cascade approach (F3). The task has been that of rewriting the Schrodinger equation and, 
more importantly, that of removing the unwanted degrees of freedom, namely, the ph and the virtual pion degrees of 
freedom, which are not present in t he ca scade-like calculation of Ref . (2^] . We have then obtained a set of equations 
for the pionic matrix density, Eqs. (7.16), in which a classical expansion can be made naturally. 

We have started by formulating the relation between the S'-matrix and cross section in the phase space approach, 
(2.13), then we have shown that the analogous relation can be written using directly the evolution operator, (2.21). 
This relation is closer to a cascade model since the simulation is set just to provide the transition probability from 
an initial point in phase space to any other point at large later times. Both in quantum mechanical calculations and 
in cascade methods, the transition probability (evolution operator) is obtained by solving the associated differential 
equation (Schrodinger equation). The rest of the paper is devoted to compute the input to be used for that equation. 
For this input to be useful it has to be written in a way that connects with cascade calculations, and we have argued 
that this suggests to ca rry ou t an expansion around the classical limit. We emphasize, however, that the validity of 
the final formulas (6.14) and (7. IE) do not rely on classical-like approximations (they are not leading order terms of 
a classical expansion). They are translations of fully quantum mechanical relations, reformulated in a phase space 
form with the help of the Wigner transformation. (Nevertheless, we recall that these formulas do not implement 
indistinguishability of the ph and this can be regarded as a classical-like feature.) 

In Sections [II and IV we have introduced the necessary formalism related to the Wigner formulation. The phase 
space (space-momentum) part of this formulation has been exploited already by various authors in the literature even 

The relevance of the time- energy form of the Wigner 
IV for the elastic channel and in Section ^ for inelastic 



in the specific subject of quantum scattering [|8[ [hi], [II 
transformation for many-body problems is noted in Section 

channels (within a simplified model). This form is used ubiquitously later in the paper. In Sections [Vl| and |VII| w e have 
introduced the necessary definitions of the matrix densities until we have pinned down the quantity, / in ( |7.2[ ) , that 
can naturally be identified with the density of particles described by a cascade-like method and we have also derived 
the evolution equations satisfied by this quantity, ( 7.15[ ) . (We note that the validity of these evolution equations holds 
regardless of the interpretation given to /.) At the same time, we have systematically studied the classical limit of 
the main formulas to verify their consistency and intuitive meaning, besides, is also within a classical expansion where 
they ca n be simple enough to be of any utility. The classical-like expansion of (7.16) has been pursued further in 
Section VIII to isolate the leading contributions to the A-body quasielastic and absorption mechanisms. 
Among the conclusions of this study we note: 

1. The separation of the pion "width" in the nucleus into a quasielastic one plus an absorptive one, is achieved by 

Here, we have already been working with the 



means of the Cutkosky rules 24 



|53j] in the diagrammatic approach, 
imaginary parts of the propagators from the very beginning and with its different analytical cuts, so that no further 
separation is needed. 

2. We have been able to give a meaning to the concept of effective A-body quasielastic and absorpt ion p robabilities. 
(However, beyond lowest orders in Ti the prescription is not unique: different choices of in Eq. ( [7.20 ) could have 
been taken.) In particular, an explicit answer is given to the problem of distinguishing three-body absorption from a 
quasielastic step followed by two-body absorption. The "genuine" three-body absorption is given by A^ in (|8~23b|), 



whereas quasielastic followed by two-body absorption is described by Q^ 1,0 ' in (8.16) and A^ 3,1 ' in ( |S.23a ). The rule 
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to obtain ^4( 3 >°) is essentially to subtract A^ 3,1 ) from the full calculation of n — > (ph) 3 obtained through a proper 
Feynman diagrammatic calculation. In practice this is done by computing the relevant pion self-energy graphs and 
then applying Cutkosky rules to pick up the imaginary part, corresponding to putting the final particles on their mass 
shell. The procedure of separation is further discussed in greater detail in ||| . 

3. Higher order effective quasielastic and absorption probabilities are quantum corrections to lowest orders and 
they are not positive definite. This is a direct consequence of unitarity (conservation of number of thrown pions 
minus absorbed pions). A weaker condition, namely, unitarity for large times only, would be enough but it would 
require quantum interference which is prohibitive in a simulation. It is important to note that in physical cases the 
genuine three body absorption is by far dominated by the collision of the pion with three nucleons which exchange 
heavier mesons, rather than by exchanging far off-shell pions. As this heavy mesons are necessarily virtual, the 
subtractions discussed here have no effect and their contribution is positive definite, very much the same as in the 
two body absorption A^ 2 ^ of our model (53). As a matter of principle, the problem of negative probabilities can be 
handled by the known method of assigning weights to the particles as they cascade, in this case a negative weight. 
Unfortunately this method introduces large statistical fluctuations. (This is th e ubiquitous negative sign problem in 



quantum simulations |54| |5g, |5(|) As noted, the concrete choice of in ( 7.15 ) cannot change the result if computed 
to all orders (although it may affect the rate of convergence of the expansion). Perhaps this ambiguity can be used 
for reducing the importance of the negative regions in the higher order quasielastic and absorptions rates. 

4. Nuclear matter and the local density prescription appear naturally in this scheme. Because the effective reaction 
rates comes about by integrating virtual degrees of freedom, which are quasi-local, the nuclear matter calculation 
results as the leading order in a semi-classical expansion. However, one of the results of the study in Ref. [Q is that 
it is important to take into account the finite range of the pion-nucleon interaction which is p-wave. This appears in 
our scheme only at higher than leading terms in Ti. 

5. Indistinguishability of particles has not been implemented in this scheme. This may or may not be more than a 
technical problem. It is certainly a quantum effect. (See ]57j ] for the use of the Wigner transformation for an optimal 
treatment of the Bose symmetry of pions.) We have chosen to remove the ph's and follow only the pions as they 
cascade through the nucleus. It should be possible to do a treatment without this removal, that is, with explicit ph's 
(of course, the virtual ph's still have to be integrated out). Perhaps this could have some incidence on the problem 
related to the missing crossed graphs. In the context transport equations, the problem of antisymmetry of nucleons 



has been successfully addressed in extensions of quantum molecular dynamics |5S|, 59 . 

6. Although with pion-nucleus scattering in mind, the scheme is more general and could be of interest for other 
problems currently dealt with by Monte Carlo simulation methods. Furthermore, the scheme is exemplified with a 
simple model, but the final equations involve Green's functions which exist in any quantum many-body theory. 
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APPENDIX A 



Let us prove Eq. ( 2.10 ) for any free Hamiltonian of the form Hq = Hq(p). From the definition of S we have 



Hn, S 



= 0, then for any p: 

e -tiH«/hg~^ e itH„/h _ g e -UH /h^ e itH /hg] ( A1 ) 

In Wigner's form, it implies 

J (2ttTi) 6 

J ^fy S (x',p';y,p)e^y- x ^ n - ttAH ^/ h (A2) 
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with AH(q.p) = Ho(p + hq) — Hq(p — hq). Upon integration over x' and b we have 

[d 2 bd 3 x'S(x\p';b + x h p) = |d 3 x'd 2 6^|^5( a; ',p';6 + y || ,p) e ^ii(y l |- a; ii)/' l -^('7,nP)/' 1 (A3) 

for large t only small values of AH can survive, and AH can be approximated by v(p) -qu with v the classical velocity. 
Then the qu integration gives = x\\ + vt. But the left-hand side does not depend on i, then it cannot depend on 

£E|| either (unless y = 0, for which case there is no scattering at all). More technically, applying f dte luJt dx\\ e - *^!!" 35 " 
on both sides of (A3) we have 



= 



{5(uj)-5(u)-AH(q p p))} J dx\\ e^lf 35 " J d 2 bd 3 x' S(x' ,p';b + x h p) (A4) 



which implies (2.10) for p ^ 0. 



APPENDIX B 



We want to prove Eqs. ( |6.5| ) and ( |6.6| ). By iterating Eq. ( p.6[ ) towards lower k, the state \E„,Ei, • • • , Ek)k must be 
given by 



\E n , Ex, . . . , Ek)k 



ifiG^\E^)ihG { ^{E k )^-F 



ihG^iEv + E k )inG%- 1) (E k _ 1 )^ 



ihG^{E„ + E k + ---+ E 2 )ihG<V(E 1 )-Fi 



K + E 1 



E k )z 



(Bl) 



and for \E)q we choose the state without ph particles, |0, E). To illustrate the method it will sufficient to prove Eq. 
(Ei) just for k = 2: 



\E n ,Ei,E2)2 



\2,E)' = 



mG^(E n )mG^(E 2 )^ 

ih&P {E„ + E 2 )ihG ( ^ (£0^ 



\E„ + E!+ E 2 ) 



dE-K dE x dE 2 
2nh 2nh 2?rft 



27rh6(E - E n - Ex - E 2 )\E^,E 1 ,E 2 ) 2 



Using the analytic properties of the propagator, we find 

"zftG« o iftGgl (£) = z?iG« (25 - ffg) , k>t, 



where is the free ph Hamiltonian in the ^-th ph subspace. And similarly 



[inG^oinG^]inG^\oinGm(E) 



nG$\E-H§-H$)mGV{E-H§\ 

ihGi 2) o ih&n o ih&$\ \ih&w o ihG<v 



(E) 



(B2) 



(B3) 



(B4) 



(B5) 



where we have used the fact that do not contain operators related to the fc'-th ph if k' < k. This is due to our 
assumption that only direct graphs are included in the propagators. The same steps in (B5) go through when the 
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operators are in place, as in Eq. QB3|), because the last operator and G p ^ commute: our assumption is that each 
time F' creates the latter ph, in this case that labeled with (2). Then 



\2,E)' 



using Eq. (pJJ), 



and finally using Eq. ( ) 



in agreement with Eq. (16 



ih 



-Ft 



ihG^oihG^ 



-F^\0,E) 



\2,E)' = G 2 (E)F^G 1 (E)F^\0,E) 
\2,E)' = G 2 (E)Fl\l,E) = \2,E). 



(B6) 

(B7) 
(B8) 



APPENDIX C 



In order to prove Eq. ( 6.11 ) let us apply the defining operator in Eq. (xE) on the recurrence (3.E) 



P W K,£^) = h 2 Gi k) (E n + iu; n ) 



x / ^e~^/ h G ph (E k + Jwtjft P (fe - 1} K + Wfc,^ + E k ;E) 
xFGl h (E k - ^Lu k )Gi k) \E v - . 
It is convenient to do a Fourier transformation of in order to make explicit the dependence on uj k : 



dE k 



tr ph lim 



2ttH t fc ->+oo 



*(*) 



x , ^| tr lim / ^ e --*W^' e ^+"*)4A 

27T?l tfc^+oo / 27r7i 



The structure of the w k integral is as follows 



W 



2irh 



-^ h h 2 G ph (E + \lo) A GUE - | w ) 



(CI) 



<G ph (E k + ±u> k )F* p( k -V(t'„,E„ + E k ;E)FGl h (E k - ±u k )&£»(E* - . (C2) 



(C3) 



where u = cv k , t = t k — t' n> E = E k and .A do not commute with G p h- By using the time representation of the 
propagators 



iHG ph (E) = / dte(t)e lt{E -"" h+l ^ /h 



(C4) 



and integrating w Eq. QC3) becomes 

= 

= fdre(t- - ) e <^+0(- E -ffph)/ ?l ie i (^-*) (iJ -- &ph)/ ' i 



and using the cyclic property of the trace 



lim tr ph W-" = tr ph 27r7i S(E - H ph )A 

t — > + oo 



(C5) 



(C6) 



from which Eq. (6.11) follows 
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Let us now prove Eq. ( 6.1 2| ) or equivalently, 

/ ^/'K,E,;E)=0. 



(C7) 



Our starting point is the Schrodinger equation ( 6.1 3| ) . There the part with ij£ , 7?^^ contains the annihilation of 
the pion from p^ k \ which go to p( k+1 \ while the other part contains the transition k — 1 — > k. Then it is enough to 
prove that the pions that disappear in pW appear in p^ k+1 ^: Eq. (|C7|) is a consequence of 



2nh 



tr 



Hl k) (E„ + \^)p^{u^E^-E)-p^\^,E v -,E)H^{E v - 



Ft E' 7r ;E)F G (k+1)t (E„ - 



= 0. 



By using the cyclic property of the trace and exchanging E n , E' n in the second part, Eq. ( |Cq ) follows from 



which is equivalent to Eq. (B4). 



(C8) 



(C9) 



APPENDIX D 



Here we wish to write the corresponding equations when an absorption mechanism is included in the model of Eq. 



( p.2|) . We simply add a 7r-ph "vertex" in Hi, 

H Trph = K + £t = J d 3 xd 3 yK(x,y)4>l(x)4> ph (y)+h.. 



(Dl) 



The vertices K , K' indicate that a pion can transform into a ph and vice versa. Combined with quasielastic, this 
implies that even if we start with one pion, time evolution will produce states with no pions, and states with many 
pions. Below the threshold this many-pion (more than one pion) stat es c an only be virtual and the exposition is 
greatly simplified by not considering them. The equations generalizing fl5.6|) are 



(D2) 



E\n(ph) k ) = H \n(ph) k ) + F^iriph)^ 1 ) + F\7r(ph) k+1 ) + K\(ph) k+1 ) 
E\(ph) k ) = ^ |(ph) fc }+Xt| 7r (ph) fc - 1 ). 



By following the same steps as before, and again keeping only direct graphs, it can be seen that all the formulae 
related to states |7r(ph) fe ) remain unchanged, except (6.4t): 



ff« = ihF \6l k+ V o G p k h +1) ] Ft + KG [k + 1] K^ 



(D3) 



the pionic Hamiltonian now contains ph's as self-energy. In this sense, the |7r(ph) fc ) states are autonomous, however, 
they act as a source for |(ph) fe+1 ) states: 



p (k \tu k ,E k ;E) = G ph {E k + \w k )fr p^ (cj k , E k ; E)K &L(E k - \uj k ) 



(k) 

where p A contains the degrees of freedom of the fc-th ph only (recall that the other k — 1 ph's are traced in p 
which is purely pionic, and we include direct graphs only). As usual, we can write an evolution equation for pA, by 
computing G~^pa ~ PaG^ 1 . Again taking trace over ph and integrating out E k , we obtain (using the trace cyclic 
property) 

ihd t N A k) {t,E) = tr e -*-*(*-*')/« 
J znTi 2-kTi 

xK^p( k ~V(t',E k ;E)K (Gl h (E k - \uj k ) - G ph (E k + \u> k )) (D5) 



(D4) 



s(fc-l) 
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where N^(t,E) = J dE k /2nh tr ffi (t, E k ; E) denotes the number of pio ns a bsorbed by k ph until time t. The total 
number of pions absorbed is obtained by integrating over t both sides of (D5) 



ihNf(oo,E)= J dttv J ^-fr p^ 1 \t,E k :E)k2Tn8{E k - H ph ) 
from which we read off the effective absorption rate 



ihd t N ( A k) 



(t,E)=ix J p {k - 1 Ht,E k ;E)k2niS(E k -H ph ). 



(D6) 



(D7) 



The difference between the exact, (D5), and the effective, (D7), absorption rates are due to quantum fluctuations 
which do not contribute to the final cross section. The simplification occurs because the ph's are not allowed to 
further interact after their creation in a quasielastic step, and so they can directly be taken on their mass-shell (as 
comes out of the formula). For pions such a simple result does not follow because the pions may always have further 
quasielastic (or absorption) steps. 

Note the odd counting in powers of fi in Eq. (D7). However, the limit fi — > will be meaningful assuming that K 
is of order h. The reason for this is that actually the ph state is not elementary, as in our model, but rather it is 
formed by a nucleon and its hole, and as a consequence the field (f>ph(x) is composite containing one loop which gives 
a Ti factor to K. 



The equation similar to (7.9b) will be 



d t N A = Pp 



where the superoperator P is given in (D7), and the equation similar to (7.15) for absorption is 

8 t N A = P (1 - PvMY 1 af Af . 



(D8) 



(D9) 



Paralleling the case of quasielastic, this formula indicates that after integration of virtual pions there will be effective 
iV-body absorption mechanisms A^ k ' k > for the real pions. 
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Assume, for instance, that p(x,p), constructed out of some p, is non-zero only in two localized regions of the phase space 
such that in one of them it is positive and is negative in the other. The positive region should be larger than C((27rfi) 3 ) 
due to the uncertainty principle while the negative one can be at most of size C((27r?i) 3 ) since p is a positive operator. If 
we were to reconstruct p using a different value for % it could not be very different in order to satisfy both constraints. 
To see this, consider arbitrarily changing the phases of \tf},t) for all ti mes t before some to- Clearly this manipulation does 
not affect p(t) for t < to, however, due to the non-locality in time in (4.3), the evolution of \ip,t) and p(t) will be modified 
after to. This shows that p(t) does not satisfy an autonomous evolution equation. 

Another simplification has been done in Appendix ^ for absorption, this time just for obtaining a simpler presentation, 
namely, to disregard ph self-energy graphs. 

As a technical remark, note that for fairly general test function spaces, the principal value of 1/x is well defined as a 
distribution, thus, unlike P(l/x k ), the construction P(l/x k ) is also a well defined distribution for k = 2,3, . . . (finite when 
applied to test functions) since the derivative of a distribution is again a distribution, defined through by parts integration 
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[66] Using E-n ± ^lo-k as arguments in the pion propagator in Eq. (7.20) (in frequency representation), would introduce further 
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terms in the definition of a^ k \ (as in ( 3.10) compared to ( 3.14b)). With this alternative definition p^) (t, E) would no longer 
be proportional to 8(E) and, as a consequence, A^ 2 ' 1 ' would not vanish and A^ 2 ' ' would not be positive at leading order 
(their sum, however, would no change). 

This is the case when Cutkosky rules [ p7| are directly applied in nuclear matter. 



